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Abstract 



The directed and elliptic flow of protons and negative pions in 4.2A GeV/c 
C+C collisions is studied using the Fourier analysis of azimuthal distributions. 
It is found that the protons exhibit pronounced directed flow, while the flow of 
pions is either non existent or too weak to be detected experimentally. Also, 
it is found that in the entire rapidity interval the elliptic flow is very small 
if not zero. These results are confirmed by the Quark-Gluon-String Model 
(QGSM) and the relativistic transport model (ART 1.0), except that these 
models predict very weak antiflow of pions. The more detailed comparison 
with the QGSM suggests that the decay of resonances and rescattering of 
secondaries dominantly determine the proton and negative pion flow at this 
energy. 

PACS number: 25.75.Ld 
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Event anisotropy, often called flow, has been observed in heavy-ion collisions at all inci- 
dent energies At Bevalac energies and below, the flow is usually studied in terms of 
the mean in-plane component of transverse momentum at a given rapidity, {p^{y)), 0] and 
additionally quantifled in terms of derivative at the midrapidity Fy = d{p^)/dy. At high 
energies, the Fourier expansion of the azimuthal distribution of particles is used [§-|10|. In 



this expansion the flrst harmonic, Vi, quantifles the directed flow while the second harmonic, 
V2, quantifles the elliptic flow. Using the Fourier expansion, the anisotropic transverse flow 
was analysed for heavy symmetric systems at the AGS 0] and SPS energies. It was 
found that this anisotropy, and particularly the elliptic flow, plays an important role for in- 
vestigating properties of hadronic matter [pJ^-p^ . However, it is still not clear whether the 
experimentally observed event anisotropy is of a dynamic origin or is due to the shadowing 
of spectator matter, passing time, etc. 

In this paper the directed and elliptic flow of protons and negative pions in 4.2A GeV/c 
C+C collisions is studied using the Fourier analysis of azimuthal distributions. The analy- 
sis is performed using 9500 C+C semicentral and central collisions obtained with the 2-m 
propane bubble chamber, exposed at JINR, Dubna synchrophasotron. The data for semicen- 
tral and central collisions roughly correspond to the upper 50% of the inelastic cross-section. 
Additionally, the same type of analysis is performed using the 400000 events generated by 
the QGSM [p!5|-T7| and 200000 events generated by the relativistic transport model (ART 



1.0) 0]. For these events the same centrality criterion is applied as in experiment, leading 
to the average impact parameter ^ 2.6. In the 2-m propane bubble chamber practically all 
charged reaction products are detected. Negative particles, except identifled electrons are 
considered to be tt". Among them remains admixture of unidentifled fast electrons (< 5%). 
All positive particles with momenta less than 0.5 GeV/c are classifled either as protons or Tr"*" 
mesons according to their ionisation density and range. Positive particles above 0.5 GeV/c 
are taken to be protons, and because of this, the admixture of of approximately 18% is 
subtracted statistically using the tt"'" and vr" momentum distributions. From the resulting 
number of protons, the projectile spectators (protons with momenta p > 3 GeV/c and emis- 
sion angle 6 < 4°) and target spectators (protons with momenta p < 0.3 GeV/c) are further 
subtracted. The resulting number of participant protons still contains some 4% of deuterons 
(with p > 0.48 GeV/c) which are statistically subtracted. The admixture of tritons (with 
p > 0.65 GeV/c) is not considered. The experimental data are also corrected to the loss of 
particles emitted at small angles relative to the optical axes of chamber. The aim of this 
correction is to obtain isotropic distribution in azimuthal angle and smooth distribution in 
emission angle (both measured with respect to the direction of the incoming projectile). 

The azimuthal distribution of particles may be represented with the flrst three terms of 
the corresponding Fourier expansion 

;^[l + 2i;icos(0) + 2t;2cos(20)], (1) 

where the two coefficients, vi and V2, quantify the directed and elliptic flow via vi = (cos^cp)) 
and V2 = (cos (20)). In Eq. (1), = (j)iab — ^ plane is the particle azimuthal angle determined 
with respect to the reaction plane, with (piat denoting the azimuthal angle of particle in the 
laboratory frame and ^piane denoting the azimuthal angle of the (true) reaction plane. Since 
both the projectile momentum and the impact parameter vectors are available in the QGSM 
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simulation, they are used to determine the corresponding reaction plane. In the experiment 
the reaction plane is determined, for each event, using the projectile momentum vector and 
the vector Q determined from |^ 

Q = Y.PTi{y > Vera + 5)- X]PTi(Z/ < Vcm " 8), (2) 
i 3 

where Pt represents the transverse momentum of the proton emitted in the forward {y > 
Vcm + S), or backward {y < ycm — S), hemisphere. Here, ycm denotes the center of mass 
rapidity while the quantity 6 (=0.2) removes the protons emitted around the ycm which 
are not contributing to the determination of the reaction plane. The reaction plane angle 
for a proton is determined using this expression only if this proton is not included in the 
above sum (i.e. if its rapidity lies in the interval from ycm — S to ycm + S)- Otherwise, in 
order to avoid autocorrelation (which is an effect of the finite multiplicity), the Q vector is 
constructed by the analogous expression in which the contribution of this proton is simply 
omitted 0. We found that the reaction plane angle distribution is essentially flat, thus 
confirming the absence of significant distorsions which could influence the magnitude of the 
extracted flow parameters. 

The accuracy with which the reaction plane angle is determined, i.e. the reaction plane 
resolution, is evaluated by the subevent method 0. In this method, each event is divided 
randomly into two subevents, and then the corresponding two reaction planes are deter- 
mined. Subsequently, the absolute value of the relative azimuthal angle, <I>i2, between these 
two estimated reaction planes is obtained. The width, a, of the $12 distribution determines 
the reaction plane resolution. For C+C collisions we find a = 50". The relative azimuthal 
angle distribution is the basis for the correction of the Fourier coefficients, v'^, obtained with 
the estimated reaction plane. The relationship between the v'^, and the Fourier coefficients 
Vn obtained relative to the true reaction plane, is v'^ = f„ (cos(nA$)), where (cos(nA$)) 
is the correction factor determined from $12 distribution following the prescription given in 
T0| , p!9[| . We find (cos(A$))=0.56 and (cos(2A$))=0.24. The correctness of this procedure 



is checked using the QGSM. Using this model, the coefficients f 1 and V2 are calculated with 
respect to the true reaction plane and also with respect to the estimated reaction plane. 
The results of the comparison will be discussed below. 

Fig. 1 (top) displays the experimentally determined f 1 coefficient vs. y (with y calculated 
in the center-of-mass frame) , for protons and negative pions together with the f 1 calculated 
with QGSM relative to the true reaction plane and relative to the estimated reaction plane. 
For the proper comparison with the experiment, we excluded protons satisfying cuts for the 
proton spectators in the experiment. In the case of protons it is seen that the values of the 
two QGSM results for t>i are quite close. The dependence of Vi on rapidity is characterised 
by a curve with a positive slope and with the zero-crossing at y = 0. The curve indicates 
a positive directed flow with magnitude vi ~ 0.17, at rapidities close to the beam rapidity 
(0.7 < y < 1.5). The QGSM reproduces satisfactorily the shape of Vi{y) curve and the 
magnitude of the flow. Using the extracted values of Vi and their relation to the mean 
transverse momentum projected onto the reaction plane, Vi = {Px)/{pt), we can evaluate 
{px) as a function of rapidity and determine the slope, F = d{px) / d{y / yb) , with respect to 
rapidity normalised to beam rapidity in the lab frame (|/b=2.2). In the present analysis we 
flnd for the slope at the midrapidity F = 144 MeV/c. After the normalisation to the mass 
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number of the colliding system we obtain the so called scaled flow Fs = F/{Ay^ + A2^^) = 
31 MeV/c. This value is in agreement with the observed trend that after reaching 
the maximum at beam energy around 0.7-2A GeV, the directed flow slowly decreases with 
increasing beam energy. 

For negative pions the experimental values of Vi indicate that the directed flow is non- 
existent. This result is confirmed by the model calculations of Vi with respect to the es- 
timated reaction plane. However, the model calculations of Vi with respect to the true 
reaction plane show the existence of a very weak directed flow of pions with the sign of Vi 
opposite to that of protons and with the maximum value of 0.02 around the target rapidity. 
This further suggests that in the collisions of light nuclei, like C-l-C, the very weak flow, if 
it exists, is not measurable because of the limited accuracy in determination of the reaction 
plane. 

Fig. 1 (bottom) displays the experimentally determined ^ coefficient vs. y for protons 
and negative pions. This coefficient is not corrected to the reaction plane resolution since 
the comparison of the model calculation of v'2, obtained as in the experiment, and the 
model calculation of V2, relative to the true reaction plane, indicates that the corresponding 
correction procedure for V2, as outlined above, is not applicable. The reason for this is the 
lightness of the colliding nuclei and the smallness of the elliptic flow. The uncorrected values 
of V2 show that in the entire rapidity interval the elliptic flow is small (|f2| < 0.02) if not 
zero, and this is consistent with the predictions of QGSM. The experimental values for V2 
also show that, for both protons and pions, the elliptic flow depends on rapidity and that 
around the beam and target rapidities it is positive for protons and negative for pions. The 
positive sign for protons indicates an enhanced emission in the reaction plane, while the 
negative sign for pions indicates an enhanced emission perpendicular to the reaction plane. 
This behaviour points out to the shadowing by the nuclear matter as the origin of the elliptic 
flow. 

Since the QGSM predictions are in fair agreement with the experimental results at 4.2A 
GeV/c, we use this model to clarify the question which of the processes are responsible for 
the flow effect. In this model, in collisions of light C+C nuclei, approximately 40% of protons 
and ~ 70% of 7i~ originate from decay of the lowest-lying resonances (A's, g,uj,ri and rj ). 
The rest originates from the 'non- resonant' primary and secondary interactions of the type: 
A^A^ A^A^TT, AA^ AA^, ttN irN, nNN NN. The protons and pions from primary 
interaction escape the collision zone without further rescattering and comprise ~ 5% of the 
total. Therefore, according to QGSM, we separately evaluate the flow of protons and pions 
originating from the following sources: (z) decay of resonances, {ii ) primary non-resonant 
interactions {Hi ) and secondary non-resonant interactions. 

Figure. 2 (top) shows f 1 vs. rapidity for protons and negative pions originating from 
decay of resonances, and from primary and secondary non-resonant interactions, as well 
as the overall Vi for protons and vr". (In these model calculations, the experimental cuts 
for the proton spectators were not applied, and this leads to a small difference between 
the two curves for overall Vi for protons in Figs. 1 and 2). The protons originating both 
from the decay of resonances and from the secondary interactions show the directed flow of 
similar intensity. The same applies to the antiflow of pions. The protons from the primary 
interactions show a relatively flat vi {y) distribution, while the pions from these interactions 
show a strong directed antiflow with magnitude ~ —0.13 around beam rapidity. The antiflow 
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of pions can be also explained by the shadowing effect, but the shadowing matter is different 
for pions from primary and secondary interactions since these pions are produced at different 
collision stages. Initially (at times less than the passing time, tp=4.2 fm/c) the pions are 
shadowed by the cold spectators. Later, after the spectator matter leaves the coUision zone, 
the pions are shadowed by the participant nucleons. This may be the underlying mechanism 
that leads to the different behaviour of f i for pions. In heavy nuclei collisions, additionally 
generated by the QGSM, the protons from primary interactions, similarly to the case of 
pions from primary interactions, show strong antiflow caused by the shadowing from the 
cold nuclear matter. In the coUisions of hght nuclei this shadowing is small and there is no 
preferential emission of the protons. 

Figure.2 (bottom) shows V2 vs. rapidity for protons and negative pions originating from 
decay of resonances, and from primary and secondary non-resonant interactions, as well as 
the overall Vi for protons and 7r~. It is seen that the particles from secondary interactions 
and from decay of resonaces, exhibit similar behaviour. The particles from the primary inter- 
actions show a clear negative elliptic flow, and this out-of-plane emission can be attributed 
to the shadowing by the cold spectators. 

In order to establish a less model dependent picture, the results of the experiment are 
also compared with the relativistic transport model, ART 1.0. These are shown in Fig. 3, 
where the calculations are performed both in the cascade and in the so-called 'mean field' 
mode. The cascade mode underestimates the magnitude of the proton flow (some 20%), and 
predicts a small directed antiflow of negative pions with magnitude vi ~ 0.04. In the mean 
field mode the model increases the magnitude of the proton flow. Also, the ART model 
predicts a very small proton and pion elliptic flow (|f2| < 0.01). 

In summary, the directed and elliptic flow of protons and negative pions in 4.2A GeV/c 
C-l-C collisions was examined using the Fourier analysis of azimuthal distributions of exper- 
imental events, and also by using the events generated by the QGSM and ART 1.0 model. 
It was found that the protons exhibit strong directed flow with magnitude Vi ^ 0.17 at 
rapidities close to the beam rapidity. The QGSM reproduces satisfactorily the shape of the 
Vi (y) curve and the magnitude of the flow. The ART model underestimates this magnitude 
in the cascade mode and increases this magnitude in the mean fleld mode. For negative 
pions the flow is either non existent, or too weak to be detected experimentally due to the 
limited accuracy in the determination of the reaction plane. The latter was suggested by 
the QGSM, where calculations with respect to the estimated reaction plane predicted non 
existant flow, while the calculations with respect to the true reaction plane predicted a small 
directed antiflow with magnitude Vi ^ 0.02. The predictions of the ART model are similar. 
Also, it was found that in the entire rapidity interval the elliptic flow is small (l^gl — 0.02) if 
not zero, and this is in agreement with the predictions of QGSM and ART 1.0. According to 
the QGSM, the two factors that dominantly determine the proton and negative pion flow, at 
this energy, are the decay of resonances and the rescattering of secondaries. The shadowing 
by the cold spectator matter affects only the flow of the particles produced at the early stage 
of the collision. 

The authors are grateful to members of the JINR Dubna group that participated in data 
processing, G. Skoro for making the QGSM code available to us, and I. Mendas for useful 
comments. 
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FIGURES 



FIG. 1. Rapidity dependence of vi and V2 for protons and it~ for 4.2A GeV/c C+C collisions: 

top- filled circles represent the experimental results for vi while the solid (dashed) line represents the 
QGSM calculation for vi with respect to the true (estimated) reaction plane; bottom- filled circles 
represent uncorrected experimental v'2 values (see text), while the solid (dashed) line represents 
the QGSM calculation for V2 iv'2) with respect to the true (estimated) reaction plane. 

FIG. 2. Rapidity dependence of vi and V2 for protons and 7r~ {solid line); for protons and 7r~ 
originating from decay of resonances (stars), primary non-resonant interactions (full circles), and 
secondary non-resonant interactions (open circles), for 4.2A GeV/c C-l-C collisions generated with 
the QGSM. 

FIG. 3. Experimental results for vi, V2 as a function of rapidity, compared with ART 1.0 
model calculations. 
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